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Abstract: This paper focuses on the problem of the degree sequence for a mixed random graph process 
which continuously combines the classical model and the BA model. Note that the number of step 



r»^^ , added edges for the mixed model is random and non-uniformly bounded. By developing a comparing 

^~^ , argument, phase transition on the degree distributions of the mixed model is revealed: while the pure 

l~:> . classical model possesses a exponential degree sequence, the pure BA model and the mixed model 

possess power law degree sequences. As an application of the methodology, phase transition on the 

ji^ ■ degree sequence of another mixed model with hard copying is also studied, especially, in the power law 

$^ ' 

jrt ' region, the inverse power can take any value greater than 1. 

1 Introduction and statement of the results 

Graph theory [51 [501 HH 121] is a rich research area that can be traced back to the problem on the seven 
bridges in Konigsberg considered by celebrated mathematician Eulcr in 1736. In 1950s, Hungarian 
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mathematicians Erdos and Renyi extended the graph theory into random environments and developed 

the classical theory of random graphs. In their paper [20], Erdos and Renyi define the random graph 

G„,M (ER model) which consists of . nodes and M randomly chosen edges of the aU (;) possible 

edges, and study the property of Gn,M as n — > oo, with AI = M{n) a function of n. At the time when 

Erdos and Renyi started their investigations of Gn.M, Gilbert |22| introduced a more fundamental 

random graph model G„.p as follows: Given n nodes, each of the I J distinct couples of nodes is 

(n\ 
linked with an edge with probability p. For M ^ p[ \ as ri — > oo, the models Gn,M and Gn,p 

are almost interchangeable and are subsequently called the classical random graph models in the 

literature. Clearly, the generation mechanism of the classical random graph is featured with several 

characteristics. First, the number of nodes is given a priori and keeps constant during the process 

of graph generation. Second, the edges are generated in a random manner. Finally, each edge is 

generated in an equal probability. 

On the other hand, in recent years complex networks have drawn a lot of attentions in disparate 
communities including statistical mechanics, computer networks, control theory, among others [H |4l 
IHl [HI [Sni HZ] . Various models involving random factors have been proposed and investigated. Among 
them, the model proposed by A.-L. Barabasi and R. Albert [7] (BA model) has been well received and 
can be described as follows. A graph with tiq nodes and ttlq edges is given at the beginning. Then 
the graph starts to evolve. At each time a new node with several new edges is added to the graph. 
While all these new edges are linked with the new node, the other node that links an edge of these 
new edges is selected from the existing nodes according to the principle of preferential attachment. 
Suppose that there are n nodes in the graph already, with dx^ being the degree of the i^^ node. The 
principle of preferential attachment asserts that the i*'' node is selected as the node that links one of 
the m edges with probability dxj X)i=i '^xi- It is shown that the degree distribution of the resulting 
graph obeys a power law. Different from the generation mechanism of Erdos and Renyi, for a random 
graph, the generation mechanism of BA is featured with the following characteristics. First, the size of 
the graph in terms of the number of nodes and edges is varying during the process of graph generation. 
The graph tends to evolve. Second, the added new edges are generated with unequal probabilities 
according the principle of preferential attachment. Obviously, the BA model can hardly be treated as 
an extension of the ER model. 

A natural question is how to reconcile the ER theory of random graphs and various models of 
complex networks and develop a coherent or modern theory of random graph and complex networks, 



this forms the first motivation of the present paper. As a useful step, it should be interesting to combine 
the distinct features of the two graph generation mechanisms described above and investigate various 
properties of the resulting graph. In this paper we will first introduce an evolving classical random 
graph model and then modify this classical model according to the principle of preferential attachment. 
The ER model can be easily modified in an evolving way as follows. Fix some constant /i > 0. Let's 
consider the following process which generates a sequence of simple graphs {G^ — (Vt, Et), t > 1}: 

Time-Step 1. Let Gj consists of vertices xq,xi and the edge (a;o,a;i). In general, (u,w) denotes 
the edge with endpoints u, v. 

Time-Step t > 2. We add a vertex xt to G^_i and then add random edges incident with x^: for 
any < i < i — 1, edge {xi,xt) is added independently with probability (p. A t)/t. 

The process {Gj : t > 1} defined above is called classical, for edges are added in an equal probability 
at any Time-Step, which coincides with the basic feature of ER model. 

It is easily observed that the classical model {G° : i > 1} is not appropriate for studying real 
world networks also. Actually, model {G° : t > 1} can be farther modified to the following BA model 
{Gt — (Vt, Et) : t > 1}, which fits the first motivation of us: 

Time-Step 1. Let Gi consists of vertices x^^xi and the edge (a;o,a;i). 

Time-Step t > 2. We add a vertex xt to Gt_i and then add random edges incident with xt: for 

any < i < t — 1, edge {xi, xt) is added independently with probability — ^ ^ 1; where d^^ (^ — 1) 

be the degree of Xi in Gt-i and et_i = \Et-i\. 

The second motivation for us to consider the above process {Gt : i > 0} is to model the www-typed 
real- world networks properly. We say a real- world network is of www-typed, if the following holds 

1. Excepting for all the isolated vertices (nodes), the network has only one connected component; 

2. There is no loop and multi-edge in the network; 

3. While a new vertex (node) is added, the number of added new edges (links) between it and the 
existing vertices is finite but unbounded; and 

4. Edges (links) are added in the preferential attachm,ent manner. 

Actually, to model the real world networks by random complex graphs, many new models (deferring 
from the ER model) have already been introduced. By studying complex graphs, various topological 



properties such as degree-distribution [71 [TOl [HI [21] , diameter [21 [H [HI [SU] , clustering [TT1[57], stability 
[51[51[T2] and spectral gap [Z] of these real- world networks have been presented. One of the most basic 
properties of real-world networks is the power law degree distribution, many new models have been 
introduced to explain the underlying causes for the emergence of power law degree distributions. This 
can be observed in the 'LCD model' [13]; the generalization of 'LCD model' due to Buckley and Osthus 
[lOj : 'copying' models of Kumar et al. [2^1; 'hard copying' models of Wu et al. [55]; the general models 
defined by Copper and Frieze [T3]; the growth-deletion models of Copper, Frieze and Vera Pj3j, Chung 
and Lu [T7] and Wu et al. 31J etc. The main difference between our model and those introduced in 
[lOl [T3l [T5l [THIITTI |25] and [31] is that, in our setting, the number of step added edges is random and 
non-uniformly bounded. Note that the 'hard copying' model introduced in [28] is also a model with 
non-uniformly bounded edge addition. Obviously, the model {Gt '. t > 1} seems to be a more proper 
candidate for modeling the www-typed real-world networks. 

Now, Let Dk{t) be the number of vertices with degree fc > in Gt and let Dk{t) be the expectation 
of Dk{t). Note that, in this paper, for any kind of random graph process, we always denote Dk{t) the 
number of vertices with degree fc > and Dk{t) its expectation. 

The first result of this paper is about BA model {Gt — (Vi, Et) : i > 1}, it follows as 

Theorem 1.1 For any < /i < 2, there exists positive constants Gi and G2 such that 

o Dk(t) Dk{t) -i , . 

Cifc-3 < liminf — ^^ < limsup— ^^ < Csfe'^ (1.1) 

t^oo t t^oa t 

for all k > 1. 

Remark 1.1 In this paper, the condition < fi < 2 is purely technical, and it is conjectured that our 
results hold for any fi > 0. 

By definition, excepting for the isolated vertices, Gt contains a unique connected component, we 
call it the giant component of Gt. Denote by Gt the giant component. The following is our result on 
E(|Ct|), the mean size of C*. 

Theorem 1.2 Assume that < /i < 2. Then for any small enough ly > 0, we have 

E(|Ct|) = (l-e-^)t + 0(i^). (1.2) 

Note that the hidden constant in 0(t^^^) only depends on v. 



Now, we present a mixed model which continuously combines the classical model {Gj : i > 1} and 
the above BA model {Gt : t > 1}. Fix some constants < a < 1 and /i, C > 0. Define random graph 
process {Gf — {Vt, Et) : i > 1} as follows. 

Time-Step 1. Let G" consists of vertices xq,xi and the edge {xq,xi). 

Time-Step t> 2. We add a new vertex xt to Gf_i and then 



1. with probability a, we add random edges incident with xt in the preferential attachment manner: 



for any < i < t — 1, edge (xi, Xt) is added independently with probability — ^ ^ li where 



d%.[t — 1) be the degree of Xi in G"_i; 

2. with probability 1 — a, we add random edges incident with Xt in the classical manner: for any 
< I < i — 1, edge (x^, Xt) is added independently with probability {Q A t)/t. 



It is straightforward to generalize the approach developed for Theorem 11.11 to prove the following 
corollary for {C^ : i > 1}, < a < 1: 

Corollary 1.3 For any 0<a<l, 0</i<2 and Q > 0, there exists positive constants G" and G^ 
such that 

Ctk-^ < liminf :^^ < limsup ^^^ < C^k-^ (1.3) 

t^co t t^oo t 

for all k> 1, where P ^1 + 2(1-' ^ ^ ^'^ 



an 

Remark 1.2 At any Time-Step i > C; ^^6 mean number of added new edges is S. '■= a/i + (1 — a)C 
and - — ^^ be the limit ratio of the number of the two kinds of edges in Gf . 

In the case of a = 0, we get the classical process {G^ ■ t > 1} parameterized by ^ > 0. Just as one 
expects, the model {Gj : i > 1} possesses a classical (exponential) degree sequence as 

Corollary 1.4 For random graph process {Gj : t > 1}, there exists positive constants C^ and G2 
such that 

for all fc > 0. 



^" ( itc J - ""^ t ^ 'TJ^^ -^ ^ ^^ lite' ^'-'^ 



Theorems 11.11 and Corollaries 11.31 and 11.41 exhibit a phase transition on the degree distributions 
of the mixed model {Gf '■ t > 1} while a varies from to 1. Note that phase transition on degree 



distributions of random graph process is first studied in the recent work "31] of Wu et al.. More 
precisely, [3lJ introduced a model with edge deletions and showed that, while a relevant parameter 
varies, the model exhibits power law degree distribution, a special degree distribution lying between 
power law and exponential^ and exponential degree distribution in turn. A numerical investigation to 
phase transition on degree distributions of networks can be founded in reference [33j . 

The rest of the paper is organized as follows. In Section 2, we give some useful estimates to e^, 
the number of edges in Gt- In section 3, we bound the maximum degree of vertex in Gt, and then 
prove Theorem 11.21 In Section 4, we establish the recurrence for Dk{t), then solve the recurrence by 
using a compare argument, and finally finish the proof of Thcorem ll.il In Section 5, we adopt the 
comparing argument developed in Section 4 to prove Corollaries 11.31 and II .41 In Section 6, we apply 
the comparing argument to study the phase transition on the degree sequence of a mixed model with 
hard copying. 

2 Estimates for et 

In this section we give some lemmas for e^, which will play important roles in the proofs of our main 
results. 

We first consider the increments of et- Let at = Ct+i — et and {^t : i > 1} be the natural u-flow 
generated by process {Gt '■ t>l]. Then 

Lemma 2.1 For all t > 1, we have 

E{at\^t)=f^ (2.1) 

and 

E(a^\^t)<{f^'^l)''k\ (2.2) 

for fc > 2. 

Proof. Let {pi : 1 < i < n}, n > 2, be a serial of positive numbers satisfying pi < i X]"=iP» ~ 1' 
and let {Xi, 1 < i < ?i} be the independent random variables with 

n 

Let X — V^ Xi. Clearly, to prove the lemma, it suffices to prove that 

i=l 

E{X) = ^ and E(X'=) < (/i V 1)'= x fc! V fc > 2. 



For fc = 1, it is straightforward to see that E(X'=) = E{X) = ^ < /i V 1. Assume that E(X™) < 
(/i V 1)™ X m! for some m > 1, then 



n n n 



E(X"+i)=E(^X,r+i = ^...^ J2 E(X,,...X,„X,„^J 



■i — 1 Zi— 1 im — 1 ■im+l — 1 



n n 



E ■ • • E E ^(^^1 • • • ^^™) + E ^(^^1 ■ ■ ■ ^»„JiE(x,„^,; 



n n 



^ E • • • E (^iE(^»i • • • ^■^™) + /''E(x,, . . . x,j) 



ii — l im— 1 



n n 



<(m+l)O.Vl)^... ^E(X,,...X,„) 



ii — 1 ini — 1 

= (m + l)(/i V 1)E(X") < (/i V l)"+i X (m + 1)!. 

Thus we finish the proof by induction. D 

Now, define Yt = et — iJd for i > 1, then, by the definition of Gt^ {Yt : i > 1} forms a martingale 
with respect to {^t : ^ > !}■ 

Lemma 2.2 There exists some constant ci > such that 

P(|ei-Mi| >i*/^)<ci<-3/5 (2.3) 

for all t > 1 . 

Proof: By the property of martingale, first, we have 

Ct-i \ 2 t-i t-i 

Y,iy^+l-y^)] =5IlE(K,+i-i;)2 = ^Var(aO. (2.4) 

i=l / 4=1 1=1 

Then, by Lemma |2. II 

t-i 
E(yt-yi)2 = ^Var(a,)< (2(^ V l)^ - /x^-) (^ _ ^^^ (2.5) 

1=1 

Finally, using the relation that ^{Y^) = E{Yt — Yi)^ + (1 — A*)^ and the Markov's inequality, we have 
P(|, - ,,| > ,./») < !W) < (2(.Vl)--.-)g^-l) + (!-,.)' , ^^,_,,. 

for some constant ci > 0. D 



Lemma 2.3 For any i^ >0, there exists constants C2,C3 > such that 

P{\et - fit\ > vt) < 026-'=='* (2.6) 

for all t > 1. 

Proof. By Lemma [2. 1[ for small A > 0, we have 

E(e^"' \^t) = l + >^l^ + 0{X^), 
then 

E(-gAet+i) = E (E(e^'=*+^'" I ^t)) ^ E (e^'='E(e^"' | ^t)) = (1 + A^ + 0{X'^))E{e^'''). 
This implies that 



A 



E(e/«') = (1 + Am + 0(A^))*-iE(e^^^) = ^^^^^^^^,^ exp{ln(l + A^ + 0{\'))t}. 
For given i^ > 0, take A > small enough such that 

4 := (^ + j/)A - ln(l + A/i + 0(A2)) > 0. 
Taking 4 = e^/(l + A^ + ©(A^)), we have 

P(et > (m + t^)t) < E{e^^')e-^^'+''^^' < ^e""^*. (2.7) 

Similarly, for some C2 , Cg > 0, we have 

P(et < (/^ - z^)t) < e(^-'')^*E(e-^'=') < c'^e''"'^'. (2.8) 

The lemma follows from I^H} and (P?5)) . D 

3 Bounding the degree and the proof of Theorem 11.21 

For times s and t with < s < i, i > 1, let d^^ (t) be the degree of vertex Xs in Gt- In this section, we 
will concentrate on the upper bound of dx^{t) and then prove Theorem II. 21 

We say an event happens quite surely (qs) if the probability of the complimentary set of the event 
is 0(t"^) for any K > 0. 

The following is our bounding for d^^ (t) . As noted in ^3T] , our result will depend on Lemma 12. 3[ 
the exponential inequality for ct- 



Lemma 3.1 For small v > and 1 < s < t, we have 

Proof. Let Xr = dx^ (t) for r = s, s + 1, . . . , t. Conditional on Xr — x and er , we have 






(3.1) 



(3.2) 



/ uix \ II U/X \ \ tix 
where B 1, j be the {0, l}-vahied random variable with P I i? I 1, ) = ^ ) = • 

Lemma |3. II follows immediately from (|2.6p . (|3.2p and a standard argument which can be found in 

the proof of Lemma 2.1 in ^SIj. D 

Remark 3.1 Because dx„{t) and dxi{t) are same distributed, Lemma \8.1\ imvlies that 



dxo{t)<t—{\ogtf, qs. 

Now, based on Lemma [5TT1 we prove Theorem 1 1.2 1 as follows. 

Proof of Theorem \l.S\ To prove Theorem 1 1.21 it suffices to show that 



n\Vt\Ct\)^e-n + 0{t—). (3.3) 

Denote by At the maximal degree in Gt . By Lemma 12. 3[ Lemma 13.11 and Remark 13.11 we have 

At 



< Lt^ 



e-t 



qs 



(3.4) 



where L be a constant independent of t. 

For large i, let's consider the probability P(at — 0), recall that at — Ct+i — et be the increment of 
et at Time-Step t + \. By equation (|3.4p . we have 



P(at = 0) = E(/a,=o) = IE(E(/a,=o I ^t)) 

E lnia,=0 I ^t) 



^t , 1 1 

— <Lt^^-^ 

et 



et 



+E E(/,,=o I ^t) 



At 1 1 

et 



— > Lt"^^^ 

et 



E E(/a,=o I -^t) 



A 



*<Lt^-i 



et 



0{t 



-10\ 



(3.5) 



The term E(/aj=o I ^t) can be expressed as 

W..^o I ^0 - rid - ^) = expjt log (l - ^ 

..J^±t^,o(t(^Y]\-e-^,o(^). ,3.) 



i—O \i—0 / } 



hence, 



^i) 



^<Li^-i) -e-^ + 0(f5^-i). (3.7) 



E ( E(/a,=o I ^ 

Thus, (IS3D and ^^ imply that 

P(at = 0) ^e-^ + oft^-^V (3.8) 

Now, by the definition of Gt^ we have 

t t-i 

n\Vt \ Ct\) = Y.^{d:,^{t) = 0) = J2P{a, = 0), (3.9) 

equation (j3.3|) foUows immediately fi'om (j3.8p and (j3.9p . D 

Remark 3.2 For any t > 1, we have 



\at = 0) = E(E(/a,=o I ^t)) = E ( J] ( 1 - 
M<^2:. (t) 



\s=0 



\s=0 ^ 

Furthermore, equation \3. 8\) implies that lim ¥{at = 0) = e^^. 

t — >oo 

For i/ie probability P{at — ^), using Iji3.4\) again, the same arguments as in kS.^S.'T^ imply that 



lim P(af = 1) = ^e"''. (3.11) 



t— i'OO 



4 The comparing Approach and The proof of Theorem 11.11 

In this Section, we develop a comparing approach to prove Theorem 11.11 We first follow the basic 
procedures in [TB] to establish the recurrence for Dk{t). By the definition of Gt, first of all, we have 
Dq{1) = 0, Di(l) = 2 and Dk{t) = for aU fc,< with k>t>l. 



10 



Now, put £>-i(i) = for alH > 1. For i + 1 > fc > and < > 1, we have 



2et 
Taking expectation and then using the basic inequality 



2et 



et<J2^ 



t{t + l) 



and the estimations given in Lemmas 12.21 and 12.31 (|4.ip imphes that 



Dkit + 1) = Dk{t) + 



k-lDk-i{t) kDk{t) , ^,,_i/5 



2 t 2 t 



+ or^/^) + /fe(i), 



(4.2) 



where /fe(i) — P(at = fc). Note that term 0(i* ''") is independent of k. We get the recurrence for 
Dk{t) as: 



W + i)^D,it) + '-^^-\Ml^oit-^'-) + Mtl 



t + 1 > fc > 0, t > 1; 



(4.3) 



i?o(l) = 0; i?i(l)-2; I?fc(t) = 0, fc > t > 1; D_i(t) - 0, i > 1. 
To solve the recurrence (|4.3p . we need a comparing argument. Note that the recurrence as (|4.3p 
with {/fc(i)} replaced by a serial of constants can be solved directly by the method developed in [T5] . 
fm and [3T|. Let 



Fk{t) ■.^Dk{t + l)-Dk{t) 



k-lDk-i{t) , kDk{t) 



hit). 



2 i 2 i 

Obviously, i^fe(i) is a determined (or known!) function in k and t satisfying 

\Fk{t)\ <Rt-^/^, y k>0, t> 1. 

For /c > 0, define 

[ ^fe(i), ifi>fc, [ A.(i), ifi>fc, 

Ak{t) = < and gk{t) = <^ 

Then, (|4.3p can be rewritten as 

Dk[t + 1) = -Dfe(i) H — h Ak{t) + gk{t), 



(4.4) 



2 f 



t + 1 > fc > 0, t > 1; 



(4.5) 



i?o(l) = 0; I?i(l) = 2; Dk{t) = 0, k > t>l- 2?_i(t) = 0, t > 1. 



11 



By the fact that /fe(i) = for t < fc - 2 and /^(fc - 1) = F{ak-i ^ k) < fik'^ for fc > 2, similar to 
(|4.4p . we have for some i?i > 

\Akit)\<Rit-'^/^, \fk>0,t>l. (4.6) 

In the rest of this section, we wih try to solve the recurrence (j4.5p for any given function serial 
{Ak(t)} satisfying (|4.6p . The lack of the existence of such limit as limf^oo fk{t) makes it difficult to 
solve (|4.5p directly. In fact, to solve (|4.5p by the known argument developed in [T5], [TB] and [3T], we 
not only need the existence of such limits, but also need a uniform speed faster than t^'^, e > 0, of 
the corresponding convergence. But this seems impossible (see the proof of Corollarv ll.4p . we have 
to develop a new method to study Dk{t). 

By Remark 13.21 limt^oo P{at = 0) = e^^, then, for some constant p > 0, 

F(at = 0) > p > 0, V i > 1. (4.7) 

For fc > 0, let 

{0, fc> 1, f Cfc-^ fc> 1, 

and ifik = < (4.8) 

p, fc = 0; [ e"'', fc = 0, 

with C ^{fiViy X 4!. Define 

0, fc> 1, i> 1, i ipk, t>k, 

Mt) = { and ^fe(t) = \ (4.9) 

^fc, fc = 0, t > 1; [ 0, 1 <f < fc. 

By Lemma [2. 11 equation (I3.10p and the Markov's inequality, we have 

Mt) < 9kit) < vkit), yk>o, t>i. (4.10) 

Now, with gk(t) in (|4.5p replaced by ipkit) and ipk{t) respectively, we get the following recurrences for 
Dk{t) and Dk{t): 



D,{t+l) = Dkit) + ^^J^-^^ + A,{t)+Mt), t+l>k>0, t>l; 

z t z t ^^-^^^ 

^o(l) = 0; ^i(l) = 2; ^^(t) = 0, fc > f > 1; ^_i(i) = 0, i > 1; 



12 



D,(t + l) = D,{t)+'^^^^~^^ + Ak{t)+Mt), t+l>k>0, t>l; 
i^o(l) = 0; A(l)-2; ^^(t) = 0, fc > f > 1; i)_i(f) = 0, f > 1. 

We first give the following comparing lemma to show that Dk{t) and Dk{t) are lower and upper 
bounds for Dk(t) respectively. 

Lemma 4.1 [Comparing Lemma] Assume that Dk{t) and Dk{t) he the solutions of Ii4.11\ ) and 
\4-12\ l respectively. Then 

Dkit) <Dkit) < Dkit), Vfc>-1, i>l. (4.13) 

Proof. We only prove the first inequality in (|4.13p . the situation for the second one is the same. 
Firstly, noticing that -D_i(t) = D-i{t) = for all t > 1, we have 

Doit + 1) = Doit) + Aoit) + Mt) 

and 

Do{t + 1) = Doit) + Aoit) + goit) 

for all t > 1. This, together with the fact that £'o(l) = -^0(1) = and the inequality (|4.10p . implies 

Doit) < Doit), Vt>l. (4.14) 

Secondly, by the fact that D^+iik) = Dk+iik) — ipk+iik) ~ gt+iik) — for all fc > 1, we have 

bk+iik + l)^-Dkik)+Ak+iik) (4.15) 



and 



Dk+iik + 1) = -Dkik) + Ak+iik) 



for all fc > 1. This, together with the initial condition 131(1) = -01(1) = 2, implies that 

Dk{k)^Dkik), Vfc>l. (4.16) 

Suppose we have proved that for some m > 0, 

Dkik + m) <Dkik + m), V fc > 1. (4.17) 

13 



If we can prove 



Dk{k + {m + l))<Dkik+im + l)), V fc > 1, 



(4.18) 



then we get the lemma by induction. 

By (|4.10p and (|4.17p . (I4.18P can be easily proved by induction. The details are omitted. D 



Now we begin to solve (j4.1ip and (|4.12p . We introduce two recurrences with respect to {4'k} and 
{(fk} as follows: 



dk = — 7^dk-i - -dk +tpk, k> 0, 



(4.19) 



fc — 1 - k - 

dk = — — rffe-i - -dk + (/5fe, fc > 0, 

(4.20) 
d_i = 0. 

The following Lemma show that (|4.19p and (|4.20p are good approximation to (|4.1ip and (|4.12p re- 
spectively. 

Lemma 4.2 Assume that {Dk{t) : fc > —l,t > 1} (resp. {Dk{t) ■ k > —l,t > 1}) be the solution of 
recurrence (OTTTp (resp. \4-l^ )and {dk : k > —1} (resp. {dk '■ k > ~l}) be the solution of Il4-19\ l(resp. 
^4-20^ ). If dk < C/k (resp. dk < C/k) for fc > and some constant C, then there exists constant Mi 
(resp. M2) such that 



bk{t) ~ tdk < Mit^l^ {resp. Dk{t) - tdk < M2t^'^ 



for all k > —I and t > 1. 



(4.21) 



Proof of Lemma \4.2\ By using the fact that Dk(t) = (resp. Dk{t) = 0) for fc > i > 1 and the 
condition dk < C/k (resp. dk < C/k)^ it is straightforward to prove Lemma 14.21 by induction (in t). 
Note that our inductive hypothesis is 

T-tt- \Qk{t)\ < Mit^/^ for all fc > -1. {resp. Ht ■ |e)fc(t)| < Mat^/^ for all fc > -1.) 



For details, one may refer to [31] (the proof of Lemma 2.2). 
Now, we finish the proof of Theorem ll.il as follows. 



n 
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Proof of Theorem \l.l\ For any given constant number serial {0^ : k > 0}, the recurrence in k with 

the form 

fc — 1 k 

dk = — -— rffc-i - i^dk + (j)k, k>0, 



can be directly solved as: (i_i = 0, do = 0o, di = |(/)i and 

7,. L. 

Applied to {V'fe} and {^k}, the summation in the right hand side of equation (|4.22p converges as 
k —> cxD, thus, dk and dk decay as k~^. Clearly, dk and dk satisfy the requirement of Lemma 14.21 and 
for some constants Ci, C2, 

Cik-^ < dk, dk < Czfc^^ V fc > 1. (4.23) 

By Lemma im Lemma [4.21 and equation (|4.23p . we have 

^,-3^7 ,. Dkjt) ^,. . fDk^ Dki^ bkjt) 7^^,_3 
Cifc < dk ^ hni < hmmt < limsup < lim — dk < C2k 

t^oo t t^oo t t^oo t i— »oo t 

for all fc > 1. n 



5 Proofs of Corollaries 11.31 and 11.41 



In this section, we prove Corollaries 11.31 and 11.41 Because the basic approach is the same as we have 
used in the proof of Thcorem ll.il we only give out a sketch. 

For the process {Gf : t > 1}, < a < 1, denote by et the number of edges in Gf and at = et+i — Ct 
none the less. 

Sketch of the proof of CoroUarv \1.3\ For simplicity, we only deal with the special case of /i = C. 

Firstly, it is straightforward to check that Lcmmas l2.HI2.2l and l2.3l hold for et- Then the recurrence 
of Dk (i) can be derived as 

— . ^ — ,^ fa(k-l) , , \Dk-i(t) fak , , \ DkU) 

Dk{t + 1) = Dk{t) + ^ , ' + (1 - a)/x ^lii _ + (1 _ a)^ ^ ^^ ^ 



2 ' ' J t V 2 J t 

+Ak{t)+g'^{t), t + l>fc>0, t>l; (5.1) 

, :Do(1) = 0; :Di(1) = 2; Dk{t) ^ 0, k > t > 1; D_i{t) ^ 0, t > I. 
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where Ak{t) satisfying ^i^, g;;^(i) == 0, V t < fc - 1 and 



j:Bh&.)^k]^a-aw(B(t^l,-J^)^k 



2et J " V V i + 1 

^■.af^{t) + (l-a)fk{t), yt>k. (5.2) 

In the case of < a < 1, we have 

k 

liminfg^(t)>(l-a) lim^(i) = (l-a)^e-^ V fc > 0, (5.3) 

then, there exists some p > such that (|4.7p holds. Note that here we get such p from (|5.3p . but in 
case of a = 1, we get it from the existence of Hmt_>oo ]P(at — 0), which depends on the degree bounds 
given in Lemma |3. II 

In case of a > 0, let n{a) =3+ [2/aJ , where [2/aJ be the integer part of 2/ a. It is straightforward 
to check that 

g^{t) < e-^ V t > 1 and g°(t) < ^-^ ^^„(^) ^ '^ , Vk>l,t>l. (5.4) 

Define {ipk} and {(pk} as 

' 0, A: > 1 




V'fe 

p, k = 0:, I e-^, fc = 0, 



with C{a) = (m V 1)"(") X n(a)!. Then define 



^fc(t) 



and 




^Pk, t>[^ + {l-a)p]Vk, 



5fcW, l<i< (^ + (l-a)M) Vfc; 



i > ( ^ + (1 - a)Ai ) V fc, 



5^(i), l<t< (^ + (l-a)^) Vfc. 



Thus we have 



Mt) < 9kit) < fkit), V fc > 0, i > 1. (5.5) 
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Let Dk{t) and Dk{t) be the solutions of the recurrences obtained from (|5.ip with g'^{t) substituted 
by ipkit) and ipk{t) respectively. Namely 



Dk{t+1) = Dkit) + ('^^i^ + (1 -«)M^ ^-^(^) 



2 ■ ^ ''J t 

+Ak{t)+ijk{t), t+l>k>0, t>l; 






^o(l) = 0; ^i(l) = 2; ^^(t) = 0, A: > t > 1; D^i{t) = 0, t > 1; 



and 



D,{t+l) = D,it) + ('^^i^ + (1 -a)M^ ""'-'^'^ 



2 ■ ^ "J t 

+Ak{t) + ipk{t), t+l>k>0, t>l; 



ak \ Dkit) 



Do{l)^0; Di{l) = 2; Dk{t) = 0, k > t > I; D^i{t) ^ 0, t > 1. 
Then Lemma l4. II holds and we have 

Dk{t) < Dk{t) < Dk{t), V fc > -1, t>l. 

Define the two recurrences with respect to {ipk} f^nd {(pk} respectively as 



and 



[ rf_i=0; 

a(fc-l) 



+ (1 - a)^l dfe_i - I — + (1 - a)^ ) dfc + </?*:, fc > 0, 



(5.6) 



^ rf_i=0. 
Then Lemma 14.21 holds, namely, under the condition that dk < C/k (resp. dk < C/k) for some 
constant C and A; > 1, there exists constant Mi (resp. Af2) such that 

<Af2t4/5^ 



i)fe(t) - tdk < Afit^/^ (resp. £>*:(*) - tdk 



(5.7) 



for all fc > — 1 and t > 1. 

Finally, it suffices to solve the recurrence in k with the form 

/a(fc — 1) \ f ak 

dk = ( ^ — - + (1 - ajfij dk-i -{^ + {l-a)fi]dk + <t>k. k > 0, 

d-i = 0, 



(5.8) 
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where {(j)k : fc > 0} be a serial of nonnegative numbers. Clearly, recurrence (|5.8p can be solved as: 

2 
d-i = 0, (io = -, — (f'Q and 
ba 






/? 



j+b 



E' 



1 



5a 



for /c > 1, 



(5.9) 



where /3 = 1 + 2/a and 6 = 2/a + 2(1 — a)fi/a. Applying to {ipk} and {(pk}, the summation term 
in the right side of equation (j5.9p converges as fc — > oo, this implies that dk, dk decay as k^^ . In 
particular, for some positive constants Cf and Cf", 

C^k-I^ < dk, dk < C^k-f^, V fc > 1. (5.10) 

Corollary [L3l follows immediately from ((5?6)) . ((57)) and ((5?T0)) . D 

Sketch of the proof of Corollary \1.4\ In the case of a = 0, the recurrence of Dk{t) can be derived 

^fc-i(i) ^Dk{t) 



as 



i?fe(i + i) = i?fe(i) + C- 



-c- 



-4fc(i)+5^(i), 



t+1 > fc > 0, t > (C-1) VI; 
:Do(1) = 0, :Di(1) = 2; Dk{t) = 0, k > t>l- D-i{t) = Q, t > 1; 
Dk>{k)=0, Dkik) = {k + 1), 0<fc'<fc, l<fc<C, 



(5.11) 



where 



Ak{t) - 



C {Dk{t) - Dk-i{t)) 



tit + 1) 

and glit) — fk{t), which is given in ()5.2|) with the parameter /i replaced by (. Note that the last line 
in (|5.1ip comes from the fact G^ is a complete graph while t < (. 

It is clear that |Afc(i)| < {2Q/t and then satisfies (14. 6|) . i.e., for some i?i > 0, 

\Ak{t)\ <Rit-^'^, Vi> 1, fc>0. 

For the term g^{t) — fk{t), we have 



lim fk{t) = f-e-^, V fc > 0; 

t-+oo fc! 



on the other hand, 



fkit - 1) = 



A ^c 



1- - 
t 



C_y t{t-i)---{t-k + i) c'' _c 



it-cy 



fc! 
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for alH > C V 2, and I < k <t, this implies that 



hit) < C(0)Vre"^, for all A; > 1 and t > 1 (5.12) 

fc! 



for some constant C(0) > 0. 

Now, by ()5.3|) . we choose p > satisfying ()4.7p and define {ipk}, {fk} as 



i'k = < 



0, k>l, f c(o)Le-^, fc>l, 

P, A: = 0; I e-^, fc ^ 0. 



Then, Corollary 11.41 follows from the comparing argument used above and the fact that 

for 4>k = Ipk and (pk respectively. D 

Remark 5.1 To get the degree distribution by the standard argument introduced in JlSf and J 161/ . 
appropriate upper bounds for At, the maximum degree, are always necessary. We point out that no 
bounds for Aj are used in our proofs of Corollaries \1.3\ and\1.4\ 



6 Application to the Hard Copying Model 

It is well known that, besides the BA mechanism, copying is another mechanism that may lead to 
power law degree sequence. The basic idea of copying conies from the fact that a new web page is 
often made by copying an old one. A kind of copying models was proposed in Kumar et al. [25j to 
explain the emergence of the degree power laws in the web graphs. These models are parameterized 
by a copy factor a € (0, 1) and a constant out-degree d > 1. At each time step, one vertex u is 
added and d out-links are generated for u as follows. First, an existing vertex p is chosen uniformly at 
random; then with probability 1 — a the i*'' out-link of p is taken to be the i'^ out-link of u, and with 
probability a a vertex is chosen from the existing vertices uniformly at random to be the destination 
of the i^^ out-link of u. It is proved in [25] that the above copying models possess a power law degree 
sequence as 4 - C'/j-(2-«)/(i-a)_ 

In this section, as an application of the comparing argument, we will introduce a new copying 
model, here we call it hard copying model. Note that another hard copying model is introduced in 
|28j , which is a mixed model of B A mechanism and hard copying mechanism. 
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For fixed < a < 1 and /i > 0, define random graph process {Gf = {Vt, Et) : i > 1} as follows 
Time-Step 1. Let G" consists of vertices xo,xi and the edge {xq,xi). 
Time-Step t > 2. 

1. with probability a, we generate vertex xt by copying an existing vertex Xi, < i < t — I from 
Vt-i uniformly at random. Note that in this case, all neighbors of xt are those of the copied 
vertex Xi; 

2. with probability I — a, we add a new vertex xt to Gf_i and then add random edges incident 
with Xt in the classical manner: for any < i < i — 1, edge {xi, Xt) is added independently with 
probability (/i A t)/t. 

As calculated in [55], for the present model, we have 

E(et)=/it + 0(t2°). (6.1) 

So et — \Et\ increase super-linearly when a > 1/2. This makes our model interesting and deferring 
from the model introduced in [25] . 

Another fact for the present model is that, in any case of a. 

At < i, V i > 1, (6.2) 

where At be the maximum degree of G" . 

In the case of a = 0, {G^ : i > 1} is just {G" : t > 1} and its degree sequence is given in 
Corollarv 11.41 (with /i in place of Q). In the case of a = 1, we get a pure hard copying model and, 
using (|6.2[) . the recurrence of Dk(t) can be derived as 

D,{t + 1) = D,{t) + (fc - 1) f %i^ - ^) , t + 1 > fc > 1; 

^ ^ ^^ (6.3) 

Z)i(l) = 2; Ufc(t) =0, fc> t > 1; Z)o(t) = 0, i > 1. 
By (|6.3p . it is straightforward to prove by induction (in t) that, there exists some Af > such that 

:Dfe(t) <M(t + l)i/^ Vfc>l (6.4) 

for alH > 1. Thus we obtain a degenerated degree distribution as follows. 
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Proposition 6.1 For any K > 1, we have 



t-oo t+1 ^ ^ 



furthermore, for any e > 0, we have 

K 



lim P I V Dk{t) > e(t + 1) I = 0. 

\fc=l / 



(6.6) 



For the case of < a < 1, using (|6.2p again, the recurrence of D}^{t) can be derived as 
Du{t + 1) = Du{t) + [a{k - 1) + (1 - a)^,] (^^^ f^) + (1 - a)/fc W, 

^ + 1 > fc > 0, t > (^- 1) VI; 

Dr^it + 1) = Dr^lt) + [a{k - I) + [l - a){t + 1)] (^^^ - f^) + (1 - «)/. W, ^^'^^ 

t+1 > fc > 0, 1 < t < (^- 1); 

:Do(1) = 0, :Di(1)-2; Dk{t) ^ Q, k> t>l- D^i{t) ^ 0, t > 1; 

where fk{t) is given in (|5.2p . 

By the comparing argument developed in Section 4, we can solve (j6.7p and obtain the following 
result. 

Theorem 6.2 For any < a < 1 and fi > 0, there exists positive constants Cf and C^ such that 

(ja^-i/a < ^j^ j^f :^M^ < lim sup ^^^^ < C^fc-i/" (6.8) 

t-^oo t f^oo t 

for all k > 1. 

Theorem 16 . 2 1 provides an interesting result: in the case of et increasing super-linearly, i.e. a > 1/2, 
the model processes power law degree sequence, furthermore, the inverse power lies in interval (1,2], 
which is never considered in previous literature. 

Remark 6.1 We note here that, except for i6.S^) . no bounds for ct and At are used in our proof of 
Theorem \6.S\ Clearly, 116. S^) holds for all models studied in this paper, and 116. 2\} implies 



Dk{t)^0, Vfc>t>l. (6.9) 

In fact, 116. 9\) is a key evidence to ensure Lemmas \4-.1\ amd \4-.2\ in the comparing argument. 
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